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ABSTRACT 



We calculate the polarization of synchrotron radiation produced at the relativistic reconfinement shocks, taking into 
account globally ordered magnetic field components, in particular toroidal and helical fields. In these shocks, toroidal 
fields produce high parallel polarization (electric vectors parallel to the projected jet axis), while chaotic fields generate 
moderate perpendicular polarization. Helical fields result in a non-axisymmetric distribution of the total and polarized 
brightness. For a diverging downstream velocity field, the Stokes parameter U does not vanish and the average polariza- 
tion is neither strictly parallel nor perpendicular. A distance at which the downstream flow is changing from diverging 
to converging can be easily identified on polarization maps as the turning point, at which polarization vectors switch, 
e.g., from clockwise to counterclockwise. 
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1. Introduction 



Magnetic fields - Polarization - Radiation mechanisms: non-thermal - Relativistic processes 



Relativistic jets are present in all radio-loud active galax- 
ies and are responsible for the bulk of their non-thermal 
emission. The spectral energy distributions of these sources 
usually show two broad components. The low-energy one, 
extending from radio wavelengths to the optical/UV band, 
and in some cases even to hard X-rays, is commonly inter- 
preted as the synchrotron radiation. This interpretation is 
supported by significant linear polarization measured rou- 
tinely in the optical (e.g., Stockman & Angel, 1978; Impey 
et al, 1991; Wills et al., 1992; Ikejiri et al., 2011), radio 
(e.g., Lister, 2001; Marscher et al., 2002; Homan et al, 2002; 
Pollack et al., 2003; Massardi et al., 2011) and other bands. 

Polarization measurements in sources dominated by rel- 
ativistic jets provide important constraints on the structure 
of the dominant emitting regions. In particular, a significant 
polarization degree indicates an anisotropic distribution of 
magnetic fields. This anisotropy may reflect a large-scale 
order in the magnetic field lines (Begelman et al., 1984) or 
may arise from chaotic magnetic fields compressed at shocks 
(Laing, 1980; Hughes et al., 1985; Jones, 1988) or sheared at 
the jet boundary layer (Laing, 1981; Matthews & Scheuer, 
1990). For purely chaotic magnetic fields compressed by 
perpendicular or conical shocks, high polarization degrees 
are found only when the polarization (electric) vectors are 
parallel to the projected jet axis ("parallel polarization"), 
while for polarization vectors perpendicular to the jet axis 
( "perpendicular polarization" ) , the polarization degrees are 
below 10% (Cawthorne & Cobb, 1990). To produce higher 
perpendicular polarization degrees, Cawthorne (2006) in- 
troduced a large-scale poloidal magnetic field component. 
Diverging conical shocks can be caused by collision of the 
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jet with a dense cloud (Lind & Blandford, 1985) or appear 
when the jet rapidly becomes overpressured with respect to 
its environment (Agudo et al., 2012). 

A different shock structure can result if a jet becomes 
underpressured. Then, the interaction between the rela- 
tivistic jet and its gaseous environment leads to the for- 
mation of so-called reconfinement or recollimation shocks 
(Sanders, 1983; Komissarov & Falle, 1997). These shocks 
deflect the jet flow and focus it on a cross section much 
smaller compared to that of a freely propagating jet. 
Nalewajko (2009, hereafter Paper I) calculated the linear 
polarization of synchrotron emission associated with re- 
confinement shocks for purely chaotic magnetic fields. He 
found that perpendicular polarization degrees can exceed 
20%. Because a reconfinement shock is an ensemble of 
conical shocks, this result appears to be in conflict with 
the work of Cawthorne &; Cobb (1990). The crucial differ- 
ence is that while Cawthorne & Cobb (1990) assumed a 
parallel upstream flow, in Paper I a spherically diverging 
upstream flow was considered. Because parallel upstream 
flows are often adopted in studies of relativistic jet polar- 
ization (Lyutikov et al., 2005; Cawthorne, 2006), we would 
like to point out that jet divergence makes a substantial 
difference in the resulting polarization. 

Reconfinement shocks with purely chaotic magnetic 
fields cannot account for the high parallel polarization often 
observed in blazars. Therefore, we extend the study pre- 
sented in Paper I to include large-scale ordered magnetic 
field components. High parallel polarization can be achieved 
by introducing a toroidal magnetic field component, but in 
general ordered magnetic fields may consist of both toroidal 
and poloidal components, forming a helical structure. In 
Section 2, we describe our model of synchrotron emission 
and polarization from relativistic reconfinement shocks, in- 
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troducing a simply parametrized family of global helical 
magnetic fields. In Section 3, we present the results demon- 
strating the effect of mixing chaotic and ordered (toroidal) 
magnetic fields (Section 3.1), and the effect of changing 
the pitch angle of the purely helical magnetic field (Section 
3.2). Our results are discussed in Section 4 and summarized 
in Section 5. Preliminary results were already presented in 
Nalewajko & Sikora (2010). 



2. Synchrotron emission and polarization 

We used the model of the structure of reconfinement shocks 
from Nalewajko & Sikora (2009) to investigate its polari- 
metric properties. Our attention is focused on the optical 
wavelengths, which in blazars are dominated by the syn- 
chrotron emission of highly relativistic electrons. The cool- 
ing time scale for these particles is very short, so that the 
emitting region is closely aligned with the sites of particle 
acceleration. 1 Relativistic shock fronts provide natural con- 
ditions for efficient particle acceleration (e.g. Achterberg 
et al., 2001). We are not concerned with details of these 
processes, it is only assumed that the relativistic electrons 
tap a constant fraction ?7 rac i of the post-shock internal 
energy. Since radiative cooling of these electrons is very 
efficient, practically their entire energy is transferred to 
the non-thermal radiation, which includes the synchrotron 
and inverse-Compton components. Although the inverse- 
Compton luminosity exceeds the synchrotron luminosity 
in radio-loud quasars, here we assume that their ratio is 
roughly constant along the jet. This is justified because 
both the magnetic energy density and the total synchrotron 
and external radiation energy densities scale roughly like 
r~ 2 (Sikora et al., 2009). Hence, a constant fraction ?7 syn 
of the energy contained in relativistic electrons goes to the 
synchrotron radiation. In this case, the power emitted in 
the form of synchrotron radiation can be determined from 
the rate of energy transfer into relativistic electrons, which 
in turn is proportional to the rate of energy dissipation. 

Let be the synchrotron radiation intensity, as 

measured in the co-moving frame of the post-shock fluid 
(denoted with a prime) . It is directly related to the emitted 
flux density, which is Lorentz-invariant, and can therefore 
be written in the external frame: 
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1 Observed frequency of the peak of synchrotron emission pro- 
duced by electrons of Lorentz factor 7 is v syu = VBj 2 4.8 MHz. 
Assuming that cooling is dominated by radiative losses due to 
the synchrotron emission, the cooling length scale, a typical dis- 
tance covered by a particle before losing a significant fraction of 
its energy, is 
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For magnetic fields strength B — 1 G and Doppler factor T> — 
10, optical-emitting electrons (f By n — 6 x 10 14 Hz) have l' cool ~ 
2 mpc. For comparison, high-frequency radio-emitting electrons 
(v syn = 43 GHz) have l' cool ~ 0.25 pc. 



where AA is the shock front area. Consider a shock front 
element covering distance Az and azimuthal angle A<j>. The 
surface area of this element is 



r s Az A(j> 
cos 6 



(2) 



where r s is the local jet radius and 9-^ s ) is the angle between 
the fluid velocity vector and the jet axis. The volume of jet 
(shocked jet) matter crossing this shock front element over 
time At is 



= P i{s) cAt sin<5 j(s) 



AA 



j(s) 



(3) 



where (3u s \ is the fluid velocity and 5u s -\ is the angle made by 
the velocity vector to the shock tangent. Its kinetic energy 
is then 



A£ k in,j(s) = (r j(s) - i)r j(s) p j(s)C 2 Av; (s ) 



(4) 



where Fj( s ) = (1 — /3j( s )) is the corresponding Lorentz 
factor and pj( s ) is the co-moving mass density. Substituting 
this into Equation 1, we can calculate $' for every shock 
front element. 

Radiation produced in the post-shock co-moving frame 
is subject to relativistic Doppler boosting into the external 
frame. Since this emitting region is fixed in the external 
frame, the transformation law for the stationary pattern 
should be used instead of a moving blob (see Appendix A 
in Sikora et al., 1997). Let k be the unit vector toward 
the observer and e the unit vector along the fluid velocity. 
Then, we have 



I(k) = ^-I'(k'). 
where 



F s (1 - (3 s k ■ e) 



(5) 



(6) 



is the Doppler factor. From the law of relativistic aberra- 
tion, we find 



k' =V s {k + [(F s - 1) (fc • e) - T S /3 S ] e} 



(7) 



Polarization of the synchrotron emission can be de- 
scribed by the polarization degree n and the positional 
angle of the electric vector xe (EVPA; "polarization an- 
gle"). The maximum polarization degree of the synchrotron 
radiation, corresponding to a uniform magnetic field, is 
n max = (3f> + 3)/(3p + 7), where p is the index of the elec- 
tron energy distribution N(j) cx 7 _p (e.g., Laing, 1980). 
When combining polarized emission from many sources, it 
is convenient to use the Stokes parameters 



Q = n cos(2 XE ) I, 
U = n sin (2 XE ) I ■ 



(8) 
(9) 



When EVPA is measured from the same reference posi- 
tional angle for all sources, the Stokes parameters are lin- 
early additive. In this work, xe is always measured from 
the projected direction of the jet axis z. To calculate the 
polarization angle, one needs to know the direction of the 
magnetic field, represented by a unit vector b. 

Following Paper I, we introduce a cartesian coordinate 
system (x,y, z), in which the jet axis is aligned with the 
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z-axis and the observer is in the xz-plane, inclined to the 
jet axis at angle 8 b s , hence 



k = [sin0 obs ,O,cos0 ob s] • 



(10) 



The post-shock fluid velocity (3 S = (3 s e is inclined to the jet 
axis at angle 6 S , hence 



e = [sin 6* s cos </>, sin 6 S sin <f), cos 6 S ] . 



(11) 



Next, we introduce an orthogonal coordinate system in the 
plane of the sky: 



K 



:,0, 



w 



k, r k„ 

1 - k'3, Y — 



The positional angle of the magnetic field vector is 



,'b'w' 
X B = arctan ( — ? 



(12) 
(13) 



(14) 



From this, we subtract the positional angle of the fluid ve- 
locity vector: 



Xc = arctan 



e • w 
e ■ v' 



The relative electric vector polarization angle is 



Xe — Xb — Xc ~2 ■ 



(15) 



(16) 



When transformed back to the external frame, the total 
flux is boosted according to Equation (5). The polarization 
degree and the polarization angle, measured with respect 
to the projection of transformation velocity vector onto the 
plane of the sky, are Lorentz-invariant. To find the absolute 
polarization angle in the external frame, one needs to add 
the positional angle of the fluid velocity with respect to the 
jet axis, which is given by 



Xc = arctan 



sin 9 S sin 6 



sin 6 S cos # b s cos <f> — cos 8 S sin o b s 



(17) 



The final formula for the electric vector polarization angle 
is 



Xe = Xb - X c + Xc - g • 



(18) 



Using a particular model of the reconfinement shock 
structure, one can calculate Stokes parameters for ev- 
ery shock front element and then integrate them to pro- 
duce simulated two-dimensional maps, longitudinal one- 
dimensional profiles, or the average signal. In the following 
sections, we consider different distributions of the magnetic 
fields and calculate the polarization of light observed at dif- 
ferent viewing angles. 



2.1. Chaotic magnetic field 

For a chaotic magnetic field of initially isotropic distribu- 
tion, considered in detail in Paper I, the symmetry is broken 
by compression at the shock front. The resulting magnetic 
field distribution can be characterized by the shock com- 
pression factor k = (pj/ps) and the unit vector normal to 
the shock front n. For a shock front element located at an 
azimuthal angle <f> and inclined to the jet axis by angle a s , 

n = [cos a s cos <j>, cos a s sin 0, — sin a s ] . (19) 

This unit vector transforms into the co-moving frame as 

n' = [cos a' s cos <fi, cos a' s sin 4>, — sin a' s ] , (20) 

where a' s is the co-moving inclination angle: 

a' s = 8 S — arctan [T s tan (9 S — a s )} . (21) 

Polarization degree of the resulting synchrotron radiation 
is given by (Hughes et al. 1985; see also Appendix A) 



n(fc') = n n 



(I-- 2 ) 


1 - (*' • n'f 
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1 - (fe' -n'f 



(22) 



The total intensity is related to the total radiation flux by 2 



(1 



1 - (k'-n'f 



2 + K 2 



$' . (23) 



The magnetic polarization vector is tangent to both the 
shock front and the plane of the sky, therefore b' cx (fe' x n'). 

2.2. Ordered magnetic field 

Models of the formation of relativistic jets involve mag- 
netic fields that are ordered on a global scale (Blandford & 
Znajek 1977; Blandford & Payne 1982; for a recent review 
see Spruit 2010; Beskin 2010). Initially of poloidal orienta- 
tion, at large distances the magnetic fields are shaped by 
the rotation of the central engine and the jet lateral expan- 
sion. Their most general form is that of a tightly wound 
helical structure, which at large distances is increasingly 
dominated by the toroidal component. Here, we introduce 
a one-parameter family of globally ordered magnetic field 
structures that can be easily implemented into the recon- 
finement shock models. 

At every point (z,r,(j)) of the jet region, the mag- 
netic field, measured in the fluid co-moving frame, has a 
toroidal component B t along the azimuthal angle <j> and 
a poloidal component along the fluid velocity vector, i.e. 
the R = Vr 2 + z 2 coordinate. The third component, which 
would be directed along the polar angle 9 = arctan(r/z), 
is neglected. The combination of these components forms a 
helix of pitch angle a-Q — arctan(£> p /i?t). 

For each magnetic field component, a relation between 
its values across the different locations needs to be estab- 
lished. First of all, the magnetic field structure is assumed 
to be axisymmetric, i.e. B p and B t are independent of <f>. 



2 The co-moving anisotropy of the total intensity has been ne- 
glected in Nalewajko (2009), but the consequences of this effect 
are not significant. 
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For a particular fluid element, conservation of the mag- 
netic flux implies the following scaling laws: B p oc R~ 2 and 
B t oc However, the dependence of the magnetic field 

components on the polar angle 9 is arbitrary. Here, a choice 
is made that for the same values of R, B p is independent of 
9 and B t oc sin 9. This scaling satisfies a condition that B t 
has to vanish at 9 = 0. The explicit formulas for the mag- 
netic field components at any point within the jet region 
are 



C t sin6> 
R 



(24) 



We relate the constants C p and C t to the value of the 
pitch angle as.m at the jet maximum radius r m at R m = 



tan as m = 



C n 



(25) 



Thus, the pitch angle at any point in the jet region is 



tan cub = (~~) tanaB, m • 



(26) 



Since r < r m , ctB,m is the minimum value of the pitch 
angle that can be achieved anywhere at the shock front. 
Eliminating C p , we calculate the magnetic field components 
before crossing the shock front 



B 



r m tan a B . 



B 



C 



(27) 



The value of Ct is not important for the resulting po- 
larization, hence the only free parameter of the magnetic 
field structure is aB, m . Across the shock front, the poloidal 
component is conserved, B pfi = _B p j, while the toroidal 
component is compressed by a factor k = (pj/p s ), i.e. 
i?t,s = B t j/n. Thus, the post-shock pitch angle is 



tanctB.s = 



tanc«B,r 



(28) 



In the coordinate system introduced in Section 2 to cal- 
culate the Stokes parameters, the magnetic field direction 
is 



hi oc B PtS e 



B t , s (p , 



(29) 



where e is the unit vector along the fluid velocity and 
4> = [— sin (j), cos </>, 0] . Since the magnetic field is locally 
uniform, the polarization degree is equal to n max and the 
total intensity is 



8" 



1 - (6' • k 



/\ 2 



(30) 



(see Appendix A). 
3. Results 

In this section, we extend the calculations presented in 
Paper I to cases involving ordered magnetic field compo- 
nents. We used the same model of the reconfinement shock 
structure as in Paper I. It is assumed that the jet plasma is 
cold and lowly magnetized, so that the jet does not accel- 
erate and expands conically, and the shock equations can 
be solved in hydrodynamical approximation. In all models 



presented here, a jet with the Lorentz factor of Tj = 10 
and the half-opening angle of 0j = 5° interacts with the 
external medium of pressure distribution p c oc r~ v with 
r\ = 0. We do not explore the dependence of the results on 
the values of these parameters, but rather on the viewing 
angle 9 b s , and on the parameters describing the magnetic 
field structure. 



3.1. Toroidal magnetic field 

For ctB,m = 0, in our model of the global magnetic field 
structure, the poloidal component vanishes everywhere, 
leaving the purely toroidal field. In this section, we inves- 
tigate the polarimetric properties of reconfinement shocks 
filled with a toroidal field component or a combination of 
chaotic and toroidal components. This combination is cal- 
culated by introducing a parameter / that describes what 
fractions of the total intensity arc emitted from each com- 
ponent: 



I — /^ordered + 0- /Kchaotic ■ 



(31) 



It can also be understood as the fraction of the total mag- 
netic energy density ub = B 2 / (8ir) contained in the ordered 
magnetic field component. Polarized intensities Q' and U' 
are combined in the same way. 

In Figure 1, we present emission maps of reconfinement 
shocks filled with toroidal magnetic fields. They can be di- 
rectly compared to the maps calculated for a chaotic mag- 
netic field, presented in Paper I, since the shock structure 
is the same. The total flux distribution is more concen- 
trated toward the jet axis. Shock outlines at the largest 
jet radius are less pronounced for 9 b s = @j/2 and barely 
visible for 9 ^ s = 0j. The longitudinal profiles of the total 
flux calculated for toroidal magnetic fields (not shown) are 
similar to those for the chaotic magnetic field. A major dif- 
ference is that for 9 b s < ©j, the central areas between the 
two emission peaks are up to ~ 30% dimmer. These differ- 
ences arise because emission from a uniform magnetic field 
is more anisotropic than emission from a two-dimensional 
chaotic distribution. In the areas close to the shock outline 
at the maximum radius, the magnetic field lines are point- 
ing more closely toward the observer, and consequently less 
emission can be seen. 

The polarization maps are very different from those for 
chaotic magnetic fields. For 9 b s = ©j, we find a very simple 
structure of parallel polarization vectors and polarization 
degrees very close to the maximum, as is generally expected 
for toroidal magnetic fields. However, the picture is more 
complicated for a smaller viewing angle. For 9 y, s = 0j/2, in 
the vicinity of the shock outline the polarization vectors are 
perpendicular to the projected shock fronts, which is very 
similar to the case of chaotic fields. But the central regions 
of the image show strong parallel polarization. Because the 
total flux is dominated by the areas close to the projected 
jet axis, their parallel polarization is expected to dominate 
the perpendicularly polarized areas. 

In Figure 2, we show the longitudinal profiles of polar- 
ization degree for several viewing angles for a toroidal mag- 
netic field. Polarization is generally parallel to the jet axis, 
with the exception of the middle region for 9 ^ s = 0j/2, 
where the maximum polarization degree is ~ 12%. Parallel 
polarization degrees reach very high values, almost ~ 60%, 
for 9 < 20j. For 9 b s = 30°, the polarization degree reaches 
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Fig. 1. Synthetic emission maps of the reconfinement shock 
with jet Lorentz factor Tj = 10, opening angle 0j = 5° and 
external pressure index r] = 0, filled with a toroidal mag- 
netic field and seen from different viewing angles o bs- The 
jet origin is marked with an empty circle, the reconfinement 
point with a filled circle. The outline of the shock fronts is 
shown with orange lines. Gray shading indicates the inten- 
sity of synchrotron radiation, while the blue bars mark the 
electric polarization vector of length proportional to the 
polarization degree. 



only ~ 9%. The polarization degree profile is very uniform 
along the jet axis for 9 b s > Oj, while for o bs = Oj/2 it is 
qualitatively similar to the profile calculated for the chaotic 
magnetic field. 

In the lower panel of Figure 2, we show the polarization 
degree profiles for an even mix of toroidal and compressed 
magnetic fields (/ — 0.5). Since the total flux (I) profiles for 
both magnetic field distributions are similar, the combined 
polarization degree closely reflects the linear combination 
used to calculate the combined polarized flux Q. Indeed, 
the polarization profiles for / = 0.5 have properties on av- 
erage between cases / = (Paper I) and / = 1 (upper 
panel of Figure 2). For o bs = ©j/2, the profile is similar 
to that for / = 1, with slightly lower parallel polarization 
maxima and higher perpendicular polarization maximum. 
For higher viewing angles, polarization is still parallel ev- 
erywhere along the jet axis. For # b s = Oj, the polariza- 
tion degree profile is less uniform, with polarization degree 
maxima at the level of ~ 42% around the jet origin and re- 
confinement point, and much lower polarization degrees in 
the middle region. For o b s = 20j and # b s = 30°, polariza- 
tion profiles are uniform with maximum degrees of ~ 21% 
and ~ 3%, respectively. Note that the highest sensitivity of 



Fig. 2. Profiles of the polarization degree along the pro- 
jected jet axis calculated for different viewing angles o bs in 
two cases: a pure toroidal magnetic field (upper panel; same 
as models shown in Figure 1) and an even mix of toroidal 
and compressed magnetic fields (lower panel). Projected 
coordinate (z/z r ) pro j equals for the jet origin and 1 for 
the reconfinement point. 



the polarization degree on / is evident for o b s = Oj and 

#obs = 20j. 

In Figure 3, we show the average polarization degree 
as a function of the viewing angle for several values of /, 
spanning between purely toroidal magnetic fields (/ = 1) 
and purely chaotic magnetic fields (/ = 0). The polariza- 
tion degree changes monotonically with / for every viewing 
angle. For the purely toroidal magnetic field, the average 
polarization degree is always parallel. The highest polar- 
ization degree, found at o b s ~ 7° ~ 1.40j, is very close to 
the maximum value II max ~ 69% (for electron distribution 
index p = 2). Around this viewing angle, the average polar- 
ization from a reconfinement shock depends very strongly 
on the value of /. One can obtain a wide range of polariza- 
tion degrees, from ~ 16% perpendicular to ~ 70% parallel 
by adjusting the balance between ordered and chaotic mag- 
netic field components. On the other hand, observation of 
very high polarization degrees from reconfinement shocks 
puts rather strong constraints on the actual viewing angle. 

For small viewing angles, 9 b s < 3°, the perpendicu- 
lar polarization signal produced by the chaotic magnetic 
field will be canceled by the parallel polarization from the 
toroidal field for / ~ 0.1, while for large viewing angles, 
6> bs > 20°, this depolarization will occur at / ~ 0.3. 
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50° 60° 



Fig. 3. Average polarization degree as a function of the 
viewing angle # b s calculated for different combinations of 
toroidal and chaotic magnetic field components. The thick 
line on the top corresponds to a purely toroidal magnetic 
field (/ = 1), the consecutive lines are calculated with steps 
A/ = 0.1 down to / = corresponding to a purely chaotic 
magnetic field. The shock structure is the same as in the 
models shown in Figures 1 and 2. The vertical line indicates 



the jet opening angle of 0j =5°. 



3.2. Helical magnetic field 

Here, we investigate the case of non- vanishing pitch angle of 
the magnetic field lines, «B,m ^ 0. Equation (28) indicates 
that the local post-shock pitch angles at different points 
of the shock front will vary from arctan (KtanaB.m) at the 
jet maximum radius r m to 90° when r s — > 0, i.e. close to 
the jet origin and the reconfinement point. Using analytical 
approximation for the shock structure (Nalewajko & Sikora, 
2009) for the flat external pressure distribution (r) = 0), the 
pitch angle can be explicitly estimated as 



tanQB.s(z) 



KtanaB.: 



4 (z/z r ) (1 - Z/Z T ) 



(32) 



In our numerical model for r\ — 0, the shock compression 
parameter n changes from 0.25 for z <C z x to ~ 0.2 for 
z — > z v . Hence, the dependence of the downstream pitch 
angle on z is not very different from that of the upstream 
pitch angle. 

In Figure 4, we show emission maps for three values of 
«B,m = 3°, 7°, 15°. They are calculated for the same shock 
structure as the maps calculated for the toroidal magnetic 
field («B. m = 0) shown in Figure 1. The viewing angle is 
fixed at (9 bs = ©j- Our first observation is that both the 
total and polarized flux maps are not symmetric with re- 
spect to the projected jet axis. This is because the observed 
brightness depends on the orientation of the magnetic field 
lines, while the projection of the helical magnetic field dis- 
tribution is not axisymmetric, even if it has an axial sym- 
metry in the three-dimensional space. As a consequence of 
this asymmetry, the integrated Stokes parameter U will not 
vanish. The average polarization angle will not be strictly 
parallel nor perpendicular to the jet axis. 

Maps of the total flux have a brighter lower rim and a 
dimmer upper rim with increasing aiB,m- The bright spots 
are shifted toward the lower rim and become less promi- 
nent. Longitudinal profiles of the total flux (not shown) 
show a systematic increase of the flux from the middle re- 
gion relative to the peak around the reconfinement point. 
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Fig. 4. Synthetic emission maps of the reconfinement shock 
with jet Lorentz factor Tj = 10, opening angle 0j = 5° and 
external pressure index r\ = 0, filled with a helical magnetic 
field of different minimum pitch angle aB, m and seen from 
the viewing angle £? bs = ©j • Dashed orange lines mark the 
turning point for polarization angle profiles in Figure 5. See 
Figure 1 for more explanation. 



For ae.m = 15°, the bright spots become very distorted 
and the upper rim becomes visible. For all values of aB jm , 
the central region of the image contributes little to the in- 
tegrated flux. 

Polarization vectors show a dramatic structural change 
with a moderate increase of the minimum pitch angle. For 
the toroidal magnetic field (as.m = 0), the polarization is 
strong and parallel across the entire image. For aB, m = 3°, 
the emission from the lower rim shows significantly lower 
polarization degrees, at one point it is almost completely 
depolarized. Polarization vectors corresponding to the main 
bright spot around the reconfinement point appear to rotate 
clockwise (CW) by ~ 20° with respect to the jet axis, while 
those around the secondary bright point around the jet ori- 
gin appear to rotate counterclockwise (CCW) by ~ 15°. 
Polarization vectors on the upper rim show little change 
compared to those for aB,m = 0. 

For QB. m = 7°, the polarization degrees decrease in the 
central region and along the upper rim, but remain very 
high around the bright spots. The polarization vectors are 
aligned with the upper rim and perpendicular to the lower 
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Fig. 5. Profiles of the average polarization degree (upper 
panel) and the polarization angle (lower panel) along the 
projected jet axis, calculated for a helical magnetic field 
with different values of the minimum pitch angle aB, m , in- 
dicated in the legend. Both panels show exactly the same 
models, three of them are shown in Figure 4. See Figure 2 
for more explanation. 



Fig. 6. Average polarization degree (upper panel) and po- 
larization angle (lower panel) as a function of the viewing 
angle, calculated for a helical magnetic field with different 
values of the minimum pitch angle ae.m, indicated in the 
legend. The shock structure is the same as in the models 
shown in Figures 4 and 5. Solid vertical line indicates the 
jet opening angle Oj =5°. 



rim. Around the main bright spot the polarization vectors 
are rotated CW by ~ 40°, and around the secondary spot 
they are rotated CCW by - 30°. 

For ttB,m = 15°, polarization degrees are significantly 
higher along the lower rim than the upper rim. Their ori- 
entation is roughly perpendicular to the lower rim and per- 
pendicular to the projected jet axis along the upper rim. 
Polarization vectors at the projected reconfinement point 
are rotated CW by ~ 60° and those at the projected jet 
origin are rotated CCW by ~ 45°. 

Not surprisingly, increasing the pitch angle in the heli- 
cal magnetic field distribution has the effect of turning the 
polarization vectors from parallel to perpendicular with re- 
spect to the projected jet axis. At the same time, the total 
emission pattern changes from that dominated by bright 
spots at both ends of the reconfinement structure to that 
dominated by one of the rims. The average polarization de- 
gree can be expected to be very sensitive to the value of 

Q!B,m- 

In Figure 5, we show the longitudinal profiles of the po- 
larization degree for aB, m ranging from to 15°. Since the 
Stokes parameter U integrated over slices of the same z pro j 
is non-zero, the polarization degree and angle are plotted 
separately. The case of ckB.m = is shown with the dashed 
line in the upper panel of Figure 2. With increasing values 
of the pitch angle, the polarization degree decreases in the 



middle part of the projected jet. At the same time, it re- 
mains very high (~ 55 — 60%) around the positions of the 
jet origin and the reconfinement point. For aB, m ~ 7°, the 
average emission becomes depolarized at (z/z r ) pro j ~ 0.62. 
We call this point the turning point and mark it with the 
vertical dashed lines in Figures 4 and 5. For higher values 
of aB.m, the polarization degree increases in the middle sec- 
tion and slightly decreases at the extremes. For aa.m = 15°, 
the polarization degree at the turning point is at the level 
of -40%. 

For small pitch angles, the polarization angle increases 
(rotates CCW) before the turning point and decreases oth- 
erwise. For QB, m < 7° the modulus of the polarization an- 
gle does not exceed 45%, i.e. Q > and the polarization 
can be described as closer to parallel than perpendicular. 
Also, the polarization at the turning point is parallel. This 
changes for aB.m > 7°, when the polarization vectors rotate 
mostly CCW with increasing Zproj and the polarization at 
the turning point becomes perpendicular. For «B,m = 15°, 
the modulus of the polarization angle is generally higher 
than 45%, hence Q < 0. 

The nature of the turning point is that the average value 
of Stokes parameter U is always close to 0. The case of 
c*B,m ~ 7° is transitional, because the average Stokes pa- 
rameter Q changes sign there. It is also evident in the emis- 
sion maps in Figure 4. Polarization vectors in the vicinity 
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of the turning point are either parallel or perpendicular to 
the projected jet axis. In the map for aB, m = 3°, the polar- 
ization is mostly parallel or very low. For ctB, m = 15°, the 
polarization is perpendicular. And for aB, m = 7°, roughly 
half of the vectors are parallel and the other half perpen- 
dicular, with relatively uniform total flux. 

In Figure 6, we show the average polarization degree 
and angle as functions of the viewing angle for the same 
values of «B, m as in Figure 5. The average polarization 
degree initially decreases with the increasing pitch angle. 
The average polarization angle is very close to for 6> bs < 
Oj and rotated in the CW direction, by a value roughly 
proportional to aB, m , for # b s > 0j. It does not depend on 
the viewing angle for 8 ^ s > 20°. 

The case of ctB, m ~ 7° is once again of special interest. 
The polarization degree for all viewing angles is lower than 
~ 11%. The polarization angle changes sharply between 
<~ for o bs = @j to ~ —70° for 9 b s = 20j and remains 
at this level for larger viewing angles. For higher values of 
the minimum pitch angle, the polarization degree increases 
and the polarization angle is generally close to 90°, i.e. po- 
larization is perpendicular to the jet axis. For aB, m = 10°, 
polarization angle deviates from 90° for very small viewing 
angles, but the corresponding polarization degrees are very 
low. For aB m = 15°, the maximum polarization degree is 
~ 40% for (Us ~ 9°. 



4. Discussion 

Polarization surveys, combined with knowledge of the inner 
jet structure inferred from VLBI imaging, can provide in- 
formation on the orientation of magnetic fields with respect 
to the projected jet axis. Optical polarization electric vec- 
tors in compact radio sources are preferentially parallel to 
the inner jet (e.g. Rusk & Seaquist, 1985). This polarization 
alignment can be explained with transverse internal shocks. 
In the reconfinement shock scenario discussed in this work, 
strong parallel polarization can only be produced in the 
presence of ordered magnetic fields that contribute at least 
half the total magnetic energy density (see Figure 3) and 
have a minimum pitch angle of less than <~ 5° (see Figure 
6). 

Polarization measured in large-scale jets between promi- 
nent knots is usually perpendicular, both in radio (e.g. 
Bridle et al., 1994) and optical bands (e.g. Perlman et al., 
2006). It indicates the presence of a strong poloidal compo- 
nent of the magnetic field. However, in expanding jet the 
poloidal magnetic field component decays faster than the 
toroidal component. In large-scale jets, the poloidal com- 
ponent of the magnetic field can arise through the velocity 
shear at the jet boundary (Laing, 1981). Cawthorne (2006) 
studied polarization of synchrotron radiation from conical 
shocks filled with a combination of chaotic and poloidal 
magnetic fields. He calculated polarized emission maps and 
compared them with resolved polarimetric VLBI map of 
blazar 3C 380. He was able to reproduce a fan-like struc- 
ture of the polarization vectors for an equal contribution of 
chaotic and poloidal magnetic field components (/ = 0.5). 
His motivation for introducing the poloidal component was 
to obtain higher degrees of the perpendicular polarization. 
Reconfinement shocks or conical shocks with diverging up- 
stream flow can produce higher perpendicular polarization 
degrees without an ordered magnetic field component. 



Lyutikov et al. (2005) studied polarization from helical 
magnetic fields in relativistic cylindrical jets. In this case, 
the pitch angle is constant along the jet and therefore their 
results are difficult to compare directly with the results 
of Section 3.2, where the pitch angle is strongly position- 
dependent. However, the dependence of the average po- 
larization degree on the viewing angle and the minimum 
pitch angle (Figure 6) is in qualitative agreement with their 
Figure 7a. One important difference is that in their model 
the average Stokes parameter (U) vanishes. In our model, 
(U) 7^ 0, since we have defined the poloidal component of 
the magnetic field along the fluid velocity unit vector e (see 
Equation 29). Only when e is parallel to the jet axis, as is 
the case for a cylindrical jet model, (U) = 0. This pro- 
vides a physical interpretation for the turning point found 
in Figure 5. It corresponds to a distance z along the jet, 
for which the post-shock velocity inclination angle 9 S = 0. 
Since 6 S is a monotonically decreasing function of z, there 
can be only one such turning point, located closely to the 
point of maximum jet radius. 

Lyutikov et al. (2005) also studied the transverse pro- 
files of polarization degree. In their Figure 9a, b they show 
that those profiles should be symmetric with respect to 
the jet axis for o b s = 1/T and strongly asymmetric for 
$obs = l/(2r) or 6> bs = 2/r. Polarization vectors close to 
the jet boundary should always be perpendicular. Using 
the emission maps shown in Figure 4, one can evaluate 
the transverse polarization degree profiles along the lines 
indicating the turning point. These models correspond to 
^obsTj = 0.9 and ^obsTs slightly lower. The Q/I profiles are 
strongly asymmetric, with the Q/I value increasing from 
the lower rim to the upper rim. Polarization vectors are 
not always perpendicular at the jet boundary, but this may 
be due to the low resolution of the maps. Our polarization 
profiles appear to be roughly consistent with their case of 
o bs = l/(2r). This may indicate that the transverse polar- 
ization profile for cylindrical jets filled with helical magnetic 
fields is very sensitive to 6 \, s for values close to 1/r. 

On the other hand, some observed radio structures in 
jets can be understood without invoking any ordered mag- 
netic field component. An interesting recent example is the 
C80 knot in the jet of radio galaxy 3C 120, reported by 
Agudo et al. (2012). It has a bow-like shape with polariza- 
tion vectors aligned perpendicularly to its outline. It was 
well reproduced by a conical shock model with compressed 
chaotic magnetic field distribution. However, this model 
can hardly explain the polarization vectors perpendicular 
to the jet axis observed immediately downstream of the 
C80 knot, even when a poloidal field component is intro- 
duced. Additional perpendicular shock waves can explain 
a parallel polarization measured even farther downstream 
(knot C99), but are inconsistent with a perpendicular po- 
larization. We propose a simple explanation of the perpen- 
dicular polarization vectors downstream of the C80 knot 
by a gradual collimation of the conical shock into a roughly 
cylindrical structure. In the absence of large-scale magnetic 
fields, the polarization vectors will be perpendicular to the 
local shock outline (Paper I). 

Our model predicts a relatively uniform distribution 
of the total brightness, especially at large viewing angles. 
Hence, it cannot reproduce very bright and compact fea- 
tures like HST-1 in radio galaxy M87. Clearly, an additional 
dissipation mechanism is required to explain its behavior, 
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in particular the high optical/UV polarization reported by 
Perlman et al. (2011). 

5. Summary 

We calculated emission maps and longitudinal profiles of 
total flux and polarization degree, and average polarization 
degree for relativistic reconfinement shocks with different 
combinations of chaotic and ordered magnetic field distri- 
butions. 

For toroidal magnetic fields (Section 3.1), the average 
polarization is always parallel to the jet axis, and the av- 
erage polarization degree can reach the maximum value al- 
lowed for a particular electron distribution. The total flux 
emission maps for 9 b s > ©j are similar to those for the case 
of chaotic magnetic fields, for small viewing angles they lack 
emission from the rims. 

For helical magnetic fields (Section 3.2), which are char- 
acterized by the value of the minimum pitch angle aB,m 
measured at the point of maximum jet radius, emission 
maps are asymmetric with respect to the projected jet axis 
and hence the Stokes parameter U ^ 0. For cus.m ^ 7°, 
the total flux is concentrated on one side of the jet axis, 
where polarization vectors align with the shock outline and 
the polarization degrees are higher than on the other side. 
Longitudinal profiles show a characteristic turning point at 
(z/z r ) pro j <~ 0.6, where average U ~ 0, regardless of the 
pitch angle. It corresponds to the point, at which the post- 
shock velocity vectors are parallel to the jet axis, as for the 
cylindrical jet models studied by Lyutikov et al. (2005). It 
separates the initial jet section, in which polarization vec- 
tors are rotated counterclockwise to the jet axis from the 
section in which the polarization vectors are rotated clock- 
wise. The average polarization angle is |xe| < 45° (Q > 0) 
for aB.m < 7°, and the opposite is generally true (see Fig. 
6). For 

c^B.m ^ 7°, the reconfinement shock appears to be 
effectively depolarized for all observers. 
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Appendix A: Polarization and anisotropy of the 
synchrotron radiation 

In Section 2, we calculated the polarization properties of synchrotron 
radiation arising from different distributions of the magnetic field. 
The total emitted power is constrained, but the intrinsic anisotropy 
of the synchrotron emission has to be accounted for. The exact results 
depend on the index p of the electron energy distribution AT (7) oc 
7~ p . However, Laing (1980) and Hughes et al. (1985) showed that 
the formulas are greatly simplified for p = 3. 



A.l. Uniform magnetic field 

For a uniform magnetic field B, we have 



I(k) = CB 2 ± , 
Q(k) = I(k)U n 
U(k) = I(k)U 

11 



: cos(2xe), 
: sin (2xe) , 



(A.l) 
(A.2) 
(A.3) 



where k is the unit vector toward the observer, B± is the magnetic 
field component normal to the line of sight, xe is the electric vector 
polarization angle, n max = (3p + 3)/ (3p + 7) is the maximum po- 
larization degree of the synchrotron radiation, and C is a constant. 
We note that B 2 ^ = B 2 — (B ■ k) 2 and Xe = Xb — 7r/2, where xb is 
the magnetic field polarization angle. We introduce orthogonal coor- 
dinates in the plane of the sky (v, w), in which the polarization angle 
is measured. The magnetic vector polarization angle is 



tanxB 
hence 
cos (2xe) 
sin (2xe) 



B w 

B v ' 



tan 2 xb - 1 _ (B ■ w) 2 - (B ■ v) 2 
tan 2 xb + 1 
-2tanxB 



tan 2 XB + 1 



-2(B ■ v)(B ■ w) 



(A.4) 

(A.5) 
(A.6) 



General formulas for the Stokes parameters are 

I(k) = C [B 2 — (B ■ k) 2 } , (A.7) 
Q(k) = CTI max [(B ■ w) 2 - (B ■ v) 2 } , (A.8) 
U(k) = Cn max [-2(B-u)(B- id)] ■ (A.9) 

The C constant can be found from the total radiation flux: 

3$ 



$ = J dn k i(k) = y cb2 



C : 



8nB 2 



(A.10) 
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A. 2. Compressed chaotic magnetic field 

Following the calculations of Hughes et al. (1985), we consider an 
isotropic magnetic field distribution 

B = Bb = B[sin 9 cos <j>, sin0 sin 0, cos 0] . (A. 11) 

It is then compressed along unit vector n = [0, 0, 1] by factor k < 1 



B = B 



— sin 9 cos 4>, — sin 9 sin <f>, cos ( 

K ft 



(A.12) 



The average Stokes parameters are found by integration over all pos- 
sible unit vectors b. For an observer located at k = [sin a, 0, cos a] 
(and plane-of-the-sky coordinate vectors v = [cos a, 0, — sin a] and 
w = [0, 1, 0]), the results are 



I(k) = ^- J dn b [b 2 - (B ■ k) 2 j 



CB 2 

= -^2- [2- (1-K 2 )sin 2 a] , (A.13) 
Q(fc) = CTW J ^ _ w)2 _ {6 w)2 j = 

= CU ™f 2 (l-« 2 ) S in 2 a, (A.14) 

CUr 
47T 



U(k) = 4 max <Kl b [-2{B-v)(B-w)]=0. (A.15) 



Hence, the polarization degree is 

, s Q(k) (l-K 2 )sin 2 a 

n ^ = T^V =IImax ^M 2^-2 ■ (A.16) 

J(fc) 2 — (1 — k ) sin ct 

in agreement with Hughes et al. (1985). The polarization angle is 
Xe = 0, since Q(k) > and U(k) = 0. The C constant is found from 

• = / d£l h I{k) = 8 " ( 9 2 + K2) CB 2 , (A.17) 
hence 

9fC 2 <I> 

C= »*Q + W (A - 18) 
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